Renormalization of the electron-phonon coupling in the one-band Hubbard model 
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We investigate the effect of electronic correlations on the coupling of electrons to Holstein phonons 
in the one-band Hubbard model. We calculate the static electron-phonon vertex within linear re- 
sponse of Kotliar-Ruckenstein slave-bosons in the paramagnetic saddle-point approximation. Within 
this approach the on-site Coulomb interaction U strongly suppresses the coupling to Holstein 
phonons at low temperatures. Moreover the vertex function does not show particularly strong 
forward scattering. Going to larger temperatures kT ~iwe find that after an initial decrease with 
U, the electron-phonon coupling starts to increase with U, confirming a recent result of Cerruti, 
Cappelluti, and Pietronero. We show that this behavior is related to an unusual reentrant behavior 
from a phase separated to a paramagnetic state upon decreasing the temperature. 
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The relevance of phonons for high-temperature super- 
conductivity has been debated since the discovery of the 
high-T c cuprates. Recently, strong renormalization ef- 
fects of the electrons near the Fermi surface, observed 
in angle-resolved photoemission in several cuprates, have 
been at least partially ascribed to phonons Further- 
more, quantum Monte Carlo simulations of the Hubbard- 
Holstein model suggest that the electron-phonon cou- 
pling shows forward scattering and no substantial sup- 
pression at large U and small dopings l2|, similar as in 
the 1/N expansion for the t-J model [3j- On the other 
hand, it has been pointed out Q that at small dopings 
the Kotliar-Ruckenstein (K-R) slave-boson approach |5j 
might yield results quite different from the 1/N expan- 
sion. Below we study the influence of strong electronic 
correlations on the electron-phonon coupling using the K- 
R approach. The quantity of interest is the static vertex 
function T which acts as a momentum-dependent, mul- 
tiplicative renormalization factor for the bare electron- 
phonon coupling. 

We consider the one-band Hubbard model on a square 
lattice with nearest and next-nearest neighbor hopping, 
t and t', respectively, 



C/^n lT n a . (1) 



For the non-interacting system the dispersion relation is 
Sk = —2t(cos(k x )+cos(ky)) — At' cos(k x ) cos(ky) , and the 
density of states has a logarithmic van Hove singularity 
at it'. For this model we want to study the influence 
of electronic correlations on the coupling of electrons to 
an external field Vi- The bare coupling has the form 
H' = J2i a n i<jVi- Writing V t — gui, H' describes also 
the interaction of electrons and atomic displacements Ui 
with coupling constant g. The linear change in the one- 
particle Green's function G(p) due to V q is, 



SG(p) 

SV a 



G(p)T(p,q)G(p + q), 



(2) 



with the charge or electron-phonon vertex T(p, q) = 
—6G~ 1 (p)/SV q . The components of the three-dimen- 



sional vectors p and q consist of a frequency and a two- 
dimensional momentum. For the calculation of the vertex 
we use the slave-boson technique of Kotliar and Rucken- 
stein 0. The basic idea of our approach Q is to calcu- 
late linear responses by linearizing the saddle-point equa- 
tions for the perturbed system about the homogeneous 
saddle-point solution. We consider only paramagnetic 
solutions. Then there arc three slave-bosons, e, p, and d, 
describing empty, singly, and doubly occupied sites, and 
two Lagrange parameters AW and A^ 2 -* enforcing con- 
sistency between slave-fermions and slave-bosons. The 
linear response to a charge-like perturbation of a given 
wave- vector can be determined by solving the 5x5 sys- 
tem of linear equations given in Considering only 



static fields Vq we find 



r(p,q) = 1 



5\W 
6V a 



Z (e p + £p+q) 



Sz 

6Va 



(3) 



The first term in Eq. is due to the explicit depen- 
dence of G~ l on V, the remaining terms are obtained 
by taking the derivative of the self-energy with respect 
to V. T does not depend on frequencies because we as- 
sumed zero frequency in q and because the saddle-point 
self-energy is frequency-independent, z is given by the 
Kotliar-Ruckenstein choice 



(e + d)p 



P 



(4) 



In the limit U — > oo our approach reduces to method 
(II) of Ref. , and we have checked that for large U we 
recover the results given in their Fig. 1. 

While it was shown that the slave-boson linear- 
response method gives very good results for the charge 
susceptibility (see, e.g., Fig. 1 of Ref. 6] for a comparison 
with exact diagonalization), it is not clear a priori how 
well it will work for the charge vertex T(p, q). As a check 
we have calculated the static vertex for a small system 
using exact diagonalization. The result for the scatter- 
ing of an electron from a state just below, to a state just 
above the Fermi surface is shown in Fig. ^ Consider- 
ing the fact that in exact diagonalization the number of 
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FIG. 1: Electron-phonon vertex T(p, q) for the scattering of 
an electron with p = (tt/2,0) by a static phonon of wave- 
vector q = (0,7r/2) (see inset, where the points denote the 
allowed p vectors and the solid line the Fermi line) for the 
Hubbard model with t' = on a 4 x 4 lattice with periodic 
boundary conditions and 5 up and 5 down electrons. The 
full line gives the result of a Lanczos calculation (with fixed 
number of electrons), the dashed line the results of the slave- 
boson linear response calculation for the same lattice. 

particles is fixed, while the slave-boson calculations are 
performed in the grand canonical ensemble, the agree- 
ment between both methods is remarkably good. This 
indicates that the slave-boson approach should work well 
at zero temperature. 

To find out how well the slave-bosons work at finite 
temperatures we compare to the quantum Monte Carlo 
(QMC) calculations of Ref. @. Fig. 4 of that work shows 
the effective electron-phonon coupling g(p, q), as defined 
in their Eq. (7), for the Hubbard model on an 8 x 8 lattice 
with filling n = 0.88, calculated at the lowest fermionic 
Matsubara frequency, for an inverse temperature of (3 = 
2. For comparison, we show in Fig. [2] the results of our 
slave-boson calculations for the same model at u> = and 
a slightly different filling n — 0.875. Also here we find a 
remarkable agreement. In particular, we find that after 
an initial decrease the coupling for forward scattering 
(small q) starts to increase for U > 8. This seems to 
indicate that the slave-boson method also works well at 
finite temperatures. Moreover, Eq. © naturally explains 
why the QMC results for different electron momenta p, 
shown in Fig. 4 of Ref. Q, are so similar. 

After comparing the results of the slave-boson linear 
response calculations to more accurate methods, which 
are, however, limited to small systems (exact diagonal- 
ization) or finite temperatures (quantum Monte Carlo), 
we now turn to very large systems at very low tem- 
peratures. First we calculate the electron-phonon ver- 
tex for electrons on the Fermi surface, where in Eq. © 
e p and £ P +q are both replaced by the Fermi energy of 
the non-interacting system. Fig. [3] shows the vertex for 
momentum transfer q along high symmetry lines in the 
Brillouin zone for the Hubbard model at essentially zero 
temperature. The effect of next-nearest neighbor hop- 
ping is illustrated by comparing calculations for t' = 



FIG. 2: Effective electron-phonon coupling for the Hubbard 
model with t' — on an 8 x 8 lattice with filling n = 0.875 
(28 + 28 electrons) as calculated in slave-boson linear response 
at an inverse temperature of j3 — 2 and for electron momenta 
on the iso-energy-line e p — (non-interacting Fermi surface 
for half- filling) . As follows from Eq. the vertex in slave- 
boson linear response is then independent of p. Thus both 
plots in Fig. 4 of Ref. can be compared to the curves shown 
above. We note that the two calculations differ slightly in the 
chosen filling and in that the QMC calculations of Ref. 
have not been done at tv = 0. 

and t' = — 0.35 1. We find that in both cases the on- 
site Coulomb interaction strongly reduces the electron- 
phonon coupling. This is not completely unexpected as 
the charge response should be strongly suppressed by an 
on-site Coulomb interaction. It is, however, in striking 
difference to the behavior at higher temperature (Fig.|2J- 
Also, while T(q) shows a broad peak around q = 0, we do 
not find a particularly pronounced forward scattering. In 
fact, the electron-phonon vertex is often strongest close 
to q = (ir, 0). This is different from what was found 
within an l/N expansion 0]. The l/N expansion relies 
on the smallness of 1/6N, i.e., it breaks down at small 
dopings S. This can be seen from the fact that the charge- 
charge correlation function remains in leading order finite 
for 5 — > though the exact correlation function vanishes 
in this limit. The Kotliar-Ruckenstein method, on the 
other hand, reproduces this limit correctly in leading or- 
der which makes it plausible that in this case the charge 
vertex is smaller than in the 1 /N expansion, especially at 
smaller dopings [4j . Which of the two methods is more 
reliable, in particular, near optimal doping, is not clear 
and can probably only be judged by comparison with 
exact numerical methods. 

To assess the importance of the electron-phonon cou- 
pling for superconductivity we calculate the renormaliza- 
tion factor 

^FS^/pS^g^PMP^-P^P') 

- 2 /fsT^t/fsT^p) 

for the pairing channels g s (p) = 1, ffs*(p) = cos(p x ) + 
cos(p y ), g Px {p) = sm(p x ), gd w2 _ y2 (p) = cos(p x )-cos(p y ), 
and gd xy {p) = sin(pa;) sin(py). A Q is equal to the ra- 
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FIG. 3: Electron-phonon vertex T(q) for scattering electrons on the Fermi surface. Calculations are for the Hubbard model on 
a 1000 x 1000 lattice with t' = (top) and t' = — 0.35 i (bottom) at an inverse temperature (5 = 500/t. Fillings and Hubbard 
U are indicated in the plots. We find that in all cases the vertex function is strongly reduced with increasing U. For large U 
the changes become small and eventually limfy^oo T(q) is reached. 



tio A q /Aq , where A Q and X a denote the dimension- 
less electron-phonon coupling constants in the interact- 
ing and non-interacting cases, respectively. To judge the 
importance of forward scattering we also calculate the 
renormalization factor for transport, 

/FS^/FS^r(p,p'-p)|v(p)~v(p')|2 

2^ /fs ^ /fs ^| v(p) | 2 

The results are shown in Fig.0| We find that for U < 10 
the s-wave couplings decrease almost exponentially with 
U . For the special case of the Hubbard model with near- 
est neighbor hopping only (t! = 0) we have A* = A s , since 
g s * is constant on the Fermi surface. Moreover A tr ~ A s , 
reflecting that there is no pronounced forward scatter- 
ing; only for larger U does A tr become somewhat smaller 
than A s . But by then both coupling constants are already 
very small. The higher pairing channels are even weaker, 
starting from zero at U = 0, going through a maximum 
around U — 2 only to decay almost exponentially. We 
can thus conclude that within Kotliar-Ruckenstein slave- 
boson theory, restricting the system to be paramagnetic, 
the contribution of Holstein phonons to superconductiv- 
ity should be very small. 

We now come back to the surprising upturn of the 
electron-phonon vertex for U > 8 shown in Fig. and 
also found in the QMC calculations of Ref . Q ■ Calculat- 
ing r(q) at kT ~ t, we indeed find a drastically differ- 
ent behavior than for T — > 0: Instead of monotonically 



decreasing with U, the coupling starts to increase and 
develops a very strong forward scattering peak. An ex- 
ample is shown in Fig. [S] Looking at the charge response 
function x(q) in the paramagnetic phase we find that 
this behavior is a precursor of a phase-separation insta- 
bility — a divergence of x(q — 0). This has already been 
pointed out in Ref. |2|. Calculating the phase diagram, 
we find a very peculiar reentrant behavior around the 
phase separated region as shown in Fig. |SJ When cool- 
ing down the system phase separates, but at low enough 
temperature it reverts back to the paramagnetic phase. 
Since in our calculations we only allow for a paramag- 
netic phase, other phases might mask the phase separa- 
tion. Also, since slave-bosons may have problems at high 
temperatures Q , it is not clear if the Hubbard model re- 
ally shows such an reentrant behavior. Nevertheless, we 
find a qualitatively similar behavior in the limit U — > oo 
in the gauge invariant 1/N expansion (i.e., a theory with- 
out Bose condensation) . Phase separation at finite T has 
also been proposed in Refs. H,0,0|. Moreover the good 
agreement with the quantum Monte Carlo calculations of 
Ref. @ suggests that our approach might indeed capture 
the relevant physics. It would therefore be interesting to 
test the phase diagram shown in Fig.HJwith QMC: A cal- 
culation for, e.g., (3=1 and U = 4. . . 8, i.e., at temper- 
atures and values of U, where QMC has little problems, 
should show clear signs of phase separation. Of course, 
these calculations should be done at uj = as the ex- 
trapolation from finite Matsubara frequencies might be 
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FIG. 4: Renormalization constants A a for different pair- 
ing channels and A tr relevant for transport for the Hubbard 
model with t = 0. Calculations were performed for lat- 
tices of increasing size at decreasing temperatures, perform- 
ing the p integrals over the whole Brillouin zone and weight- 
ing with minus the derivative of the Fermi-Dirac distribution 
~f'{ £ v) — /3/( £ p)(l — /( £ p))i which becomes a delta func- 
tions for T = 0. Convergence of the T — > extrapolation has 
been checked by comparing A s and A s * , which only for T = 
are equal. Error-bars for the extrapolation are plotted, but 
are usually smaller than the size of the plotting symbols. 
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FIG. 5: Electron-phonon vertex T(q) for scattering on the 
Fermi surface. Calculations are for the Hubbard model on 
a 100 x 100 lattice with t' = and filling n = 0.80 at an 
inverse temperature j3 = 1. The plot should be compared to 
the uppermost panel on the right of Fig. 
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FIG. 6: Phase separation: The lines enclose the region where 
the paramagnetic slave-boson saddle-point solution is unsta- 
ble against phase separation. Calculations are for the Hub- 
bard model on a 100 x 100 lattice with t' = 0. For t' /t > 
the region of phase separation tends to increase, for t'/i < 
it tends to decrease, in particular, for large doping. 



difficult close to an instability. 

In conclusion, we have studied the influence of strong 
electronic correlations on the electron-phonon interac- 
tion for the Hubbard-Holstcin model using the Kotliar- 
Ruckenstein slave boson method. For high temperatures 
the boundaries of the phase-separated state were de- 
termined in the T — U plane for different dopings and 
the increase of the static vertex T near the boundaries 
was studied, confirming and extending recent results of 
Ref. . At low temperatures and moderate or small dop- 
ings we found that T does not exhibit pronounced forward 
scattering behavior and that T reduces dramatically the 
electron-phonon coupling. It seems that exact numerical 
calculations are necessary to judge the reliability of the 
1/N and the Kotliar-Ruckenstein approaches. 
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